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SUMMARY 


The total lift responses of wings to sinusoidal gusts and to sinus- 
oidal vertical oscillations are calculated from the response to gust 
penetration and to a sudden change in sinking velocity through use of the 
well-established reciprocal relations for unsteady flow. The cases con- 
sidered are two-dimensional wings in incompressible, subsonic compressible, 
sonic, and supersonic flow; elliptical and rectangular wings in incom- 
pressible flow; wide rectangular and delta wings in supersonic flow; and 
delta wings of vanishingly low aspect ratio in incompressible and com- 
pressible flow. For most of the cases considered, closed-form expressions 
are given and the final results are presented in the form of plots of the 
square of the modulus of the lift coefficient for wings in a sinusoidally 
oscillating gust and in the form of the real and imaginary parts of the 
lift component for wings undergoing sinusoidal sinking oscillations. A 
summary table is presented as a guide to the scope and results of this 
paper; this table contains the figure and equation numbers for the types 
of flow and plan forms considered. 


INTRODUCTION 


Two of the factors required in the harmonic analysis of airplane 
response to continuous atmospheric turbulence are the unsteady-lift func- 
tions associated with sinusoidal vertical oscillations and with sinusoidal 
gusts. The unsteady-lift functions associated with a rigid wing under- 
going sinusoidal translational oscillations have been derived in refer- 
ences 1 to 10 for two-dimensional wings in incompressible, subsonic com- 
pressible, sonic, and supersonic flow; for elliptical and rectangular 
wings in incompressible flow; for wide rectangular and delta wings in 
supersonic flow; and for very narrow delta wings in incompressible and 
compressible flow. Calculations of the unsteady-lift functions associated 
with rigid restrained wings in sinusoidal gusts seem to be nonexistent, 
with the exception of the work by Jones (ref. 6) for elliptical wings in 
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incompressible flow and by Garrick (ref. ll) and Sears (ref. 12) for 
wings in two-dimensional incompressible flow. 

The purpose of this report is twofold - to compile the unsteady- 
lift functions associated with sinusoidal sinking oscillations and to 
derive the unsteady-lift functions associated with a rigid restrained 
wing in a sinusoidal gust. These latter functions are derived herein 
from existing unsteady-lift functions for a wing penetrating a sharp- 
edged gust by means of the reciprocal relation between the function for 
a wing in a sinusoidal gust and the function for a wing penetrating a 
unit sharp-edged gust. The reciprocal relation used was of the same type 
as that reported in reference 11. 

The unsteady-lift functions associated with a rigid restrained wing 
in a sinusoidal oscillating gust are derived for two-dimensional wings 
in incompressible, subsonic compressible, sonic, and supersonic flow} 
for elliptical and rectangular wings in incompressible flow; and for wide 
rectangular and delta wings in supersonic flow. In addition, the indicial 
lift function for a wing penetrating a sharp -edged gust and the corre- 
sponding oscillatory lift function are derived for a delta wing of van- 
ishing aspect ratio in compressible flow. The functions presented in 
this paper are total lift functions which include the circulatory and 
noncirculatory components . 

In studies of the airplane response to atmospheric turbulence (see 
ref. 13, for instance), the unsteady-lift functions for a rigid wing in 
a sinusoidal gust usually appear in the form of the square of the modulus 
of lift coefficient, whereas the unsteady-lift functions for a wing under- 
going sinusoidal sinking oscillations appear in the form of the individual 
in-phase and out-of-phase (real and imaginary, respectively) components 
of lift. Therefore, on this basis, all the results in this paper are 
presented in the figures in the forms mentioned. An index to the figures 
and equations or other sources of information for the unsteady-lift func- 
tions for the types of flow and wing plan forms considered herein is pre- 
sented as a table. 


SYMBOIS 


A aspect ratio 

a velocity of sound 

b (x) spanwise coordinate of leading edge of wing, measured from root 
chord, mx 

C(k) total lift coefficient for wing oscillating harmonically in pure 
translational motion, normalized to unity by its steady-state 
value, F(k) + iG(k) 
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C(k) cir 

circulatory component of C(k) 

C L 

steady-state lift coefficient 

C Ia 

wing lift-curve slope 

c 

root chord of wing 

E 

complete elliptic integral of second kind with modulus 

i - ikf 

F(k) 

in-phase component of C(k) (real part) 

f(k) 

Fresnel integral (see eq. (56)) 


Schwarz function of order n (see eq. (62)) 

G(k) 

out-of -phase component of C(k) (imaginary part) 

ho 

amplitude of vertical velocity of wing 


Bessel function of first kind 

k 

reduced-frequency parameter, coc /2V 

k^s) 

lift coefficient for wing experiencing sudden change in sinking 
speed, normalized to unity by its steady-state lift 

ki(s) cir 

circulatory component of k]_(s) 

kgCs) 

lift coefficient for wing penetrating sharp-edged gust, 
normalized to unity by its steady-state lift 

hg,ind 

lift on rigid restrained wing penetrating sharp-edged gust 

hg,OSC 

total lift on rigid restrained wing in sinusoidal gust 

^ind 

lift on rigid wing experiencing sudden change in sinking speed 

hose 

total lift on rigid wing oscillating harmonically in pure 
translational motion 

l(x) 

lift per unit length 


NACA TN 3748 


M 

m 

Ap/<1 

q. 

s 

s 

t 

t' 

V 

Wo 


x,y,z 

Y n (k) 


Mach number 

tangent of semiapex angle of delta wing 
loading coefficient 

pv 2 

dynamic pressure, 
wing area 

nondimensional distance traveled, root semichords 
distance traveled by sound wave, at 1 

time variable 
forward velocity 

amplitude of vertical gust velocity 
coordinate axes , fixed on wing 
Bessel function of second kind 


p = J|m 2 - 1| 

5(s) unit impulse function or Dirac delta function 

p air density 

cp(k) total lift coefficient for wing immersed in harmonically 

oscillating gust, normalized to unity by its steady-state 
value 

03 circular frequency 

5S = 2M 2 k/p 2 

2(s) unit jump function 


PROCEDURE 


Since this paper deals with the lift functions C(k) and cp(k), the 
lift due to sinusoidal sinking oscillations and the lift due to sinusoidal 
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gusts, respectively, a brief description of the total lift in terms of 
these functions and the method by which they were derived is in order. 

For a rigid wing oscillating harmonically in pure translational 
motion, the total lift can be expressed as 

h «= iks 

W " c(k) (i) 

where k is the reduced-frequency parameter <nc/2V, and C(k) is a 
complex quantity F(k) + iG(k). The real part of this quantity is asso- 
ciated with the in-phase component of lift and the imaginary part with 
the out-of -phase component of lift; the total lift functions include both 
circulatory and noncirculatory effects. For a rigid restrained wing in 
a sinusoidal gust the total lift can be expressed as 

W lks 

kg,osc = ^(k) (2) 

The lift on a rigid wing experiencing a sudden acquisition of verti- 
cal velocity h 0 can be expressed as 

Lind = -qSC^ ^ k x (s) (3) 

and for a wing penetrating a sharp-edged gust of vertical velocity W 0 , 
the lift can be expressed as 


L g,ind = -OSC;^ ^ k 2 (s) (4) 

where k-^(s) and ^(s) are the indicial lift functions for a wing given 

a sudden change in sinking speed and for a wing penetrating a sharp-edged 
gust, respectively- 

The functions C(k) and k-^(s) are reciprocally related as shown 
in reference 11 by the following expressions: 


C(k) = F(k) + iG(k) = 1 + ik J [k-^s) - lj e -iks ds (5) 


t (s) = 1 + -L- f Cl k ) - 1 e ikS ( 

2«i . / ^ k 


‘dk (s > 0) (6) 
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Similarly, cp(k) and k^s) are reciprocally related as follows: 

<p(k) = 1 + ik J jjsgCs) - lje -iks ds (7) 


kp (s) = 1 + — f — e iks dk (s > 0) (8) 

^ 2*iJ_ 00 k 

In appendix A the functions k L (s) and k^s) are given as obtained 

from various references for different types of flows and plan forms. The 
functions C(k) and cp(k) were obtained by means of equations (5) and 
(7), respectively, and are also presented in appendix A. The various 
types of flows and plan forms for which these functions were derived are 
discussed more fully in the following section. 


PRESENTATION OF RESULTS 


The unsteady-lift functions k-^s), k 2 (s), C(k), and <?(k) are 

presented in appendix A and in figures 1 to 20. These functions are 
given for two-dimensional wings in incompressible flow (figs. 1 and 2), 
subsonic compressible flow (figs. 7 an( 3- sonic flow (figs. 11 and 12), 
and supersonic flow (figs. 13 and l4)j for elliptical and rectangular 
wings in incompressible flow (figs. 3 to 6)j for wide delta and rectan- 
gular wings in supersonic flow (figs. 15 to 20 )j and for delta wings of 
vanishing aspect ratio in incompressible and compressible flow (figs. 9 
and 10 ) . 

The C(k) functions, although derived by other authors for all the 
wings considered herein, were recalculated by means of equation (5) from 
existing k-^(s) functions. The functions C(k) as derived by use of 

equation (5) are in agreement with the functions derived by other authors. 

The results are given by the equations in appendix A and the figures 
which contain plots of the modulus squared for the function cp(k) (that 

is, I cp(k) 1 2 ) and the separated real and imaginary parts of the func- 
tion C(k) (that is, F(k) and G(k)). As an aid to the reader, table I 
has been prepared as an index to the equation or reference identifying the 
functions k x (s), kg(s), C(k), or cp(k), the plan form and type of flow 

for which these functions were considered, and the figures where the func- 
tions |cp(k)| 2 , F(k), and -G(k) are plotted. 
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CONCLUDING REMARKS 


The total lift responses to sinusoidal sinking oscillations C(k), 
and to sinusoidal gusts cp(k), have been calculated through use of the 
well-established reciprocal relations for unsteady flow for two-dimensional 
wings in incompressible, subsonic compressible, sonic, and supersonic flow; 
for elliptical and rectangular wings in incompressible flow; for wide 
rectangular and delta wings in supersonic flow; and for delta wings of 
vanishingly low aspect ratio in incompressible and compressible flow. For 
most of the cases considered, closed-form expressions are given and the 
final results are presented in the form of plots of the square of the 
modulus of the lift coefficients for a wing in a sinusoidal gust, and the 
in-phase and out-of -phase lift components are presented for a wing under- 
going sinusoidal sinking oscillations. 

Certain gaps still exist in the knowledge of the unsteady-lift problem. 
For instance, there seems to be little or no information available for the 
swept wing. For rectangular wings in subsonic flow, and in supersonic flow 
for which the characteristic Mach lines intersect the side edges of the 
wing, the unsteady-lift problem remains unsolved, as it is for the delta 
wing for subsonic compressible and incompressible flow. Information on 
other wings with subsonic leading edges in supersonic flow is also 
missing. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 8 , 1 956 . 
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APPENDIX A 

A COMPILATION OF THE FUNCTIONS k-^s), k 2 (s), 

C(k), AND <p(k) 

Two-Dimensional Wing in Incompressible Flow 

The functions k-j_( s ) and kgCs) have been derived by Wagner 

(ref. 14) and Von Karman and Sears (ref. 15 ), respectively. Exponential 
approximations to these functions have been given by Robert T. Jones 
(ref. 6) and are now presented, together with the C(k) and <?(k) func- 
tions as given in references 1 and 11, respectively: 

k' L ( s ) « 1.0 - 0 .l 65 e"°'° 1 < ' 5 s _ o.335 e ”°‘^ 00s + \ &( s ) (9) 

kg(s) « 1.0 - 0.236e _0 '° 58s - 0.513e"°' 561< - s - 0.171e" 2,1< - 2s (lO) 

C(k) = C(k) cir + ^ (11) 

cp(k) = |c(k) cir [j 0 (k) - iJ 1 (k)J + iJpCklje - ^ (l2) 

where C(k) is defined as the total lift function and C(k) cir repre- 
sents the circulatory component of the lift and is given in reference 1 
as 


C(k) 


-J L (k) + iY x (k) 


eir 


-[j x (k) + Y 0 (k)] + i[Y L (k) - J Q (k)] 


(13) 


and 


C(k) « |(1 + ik) 

lfWl 2 ”sk 


(k > > 1.0) 


( 14 ) 


The multiplier e" ik (eq. 12) is not included in the function <p( k ) 

(ref. ll) because in reference 11 the time origin is the instant at which 
the gust reaches the midchord position of the wing, whereas in this paper 
-tqig "time origin is the instant at which the gust reaches the leading edge 
of the wing. 
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Elliptical Wing in Incompressible Flow 


For an elliptical wing in incompressible flow, the following equa- 
tions were derived from references 6, l6, and 17: 


k l (s) A=0 =i(s) + 3 5(s) 


k-, (s) » 1.0 - 0.283e"°'^ t ' 0s + 

L A=3 


86(s) 


3 C I* E 


k, (s ) « 1.0 - 0.36le"°- 58ls + 

1 A=6 


Fft 


8 B(s) 


(15) 

(16) 


5 V 


’>■ - <£[ 


(17) 


where E = E 
kind. Also, 


1 - iX 


rtA 


is a complete elliptic integral of the second 


kgCs) = s(2 - s) (0 = s = 1.0) 


k 0 (s) = 1.0 (s > 1.0) 


A=0 
-2 (S A=0 

kp(s ) « 1.0 - 0.6 T 9e-°-558 s _ 0.227e~ 5 * 20s 

^ A=3 


(18) 

(19) 


k2(s) A= 6 « 1.0 - 0.448e _0,29 ° s - 0.272e" 0 '7 25s - 0.193e" 5 -°° s (20) 


c(k) A .o - 1 + ^ 


o(k ) A=5 


ik 


_1_ _ 0.283 \ + 8ik 

ik 0-540 + ik/ 3EC 


c(k) A=6 ~ lk (ik " 0.381 + ik) T 3%^ 


8ik 


<p( k ) A=0 . ik(. 4 + 2 - 2e 


-ik' 


ik 5 


(21) 

(22) 

(23) 

(24) 
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<P(k ) A=3 


ik( JL 0-679 0.227 

\ik “ 0.558 + ik 3.20 + ik 


(25) 


<P(k) A=6 - lk (i - 0.29^+ Ik 


0.272 


0.193 


O.725 + Ik 3.00 + ik. 


(26) 


Since the functions C(k) and cp(k) have been derived from approximate 
expressions for aspect ratios of 3 and 6, no exact asymptotic expressions 
could he obtained. However, for the vanishing-aspect-ratio case, it can 
be shown that 


C(k) « 1 + ^ 

3 


(k > > 1.0) 


(27) 


(The magnitude of the Dirac delta function in equations (l6) and (17) 

\ Q 

approaches ■— as A — >00 for 

6it \ 

for s in half mean geometric chords.) 


s in half root chords but reduces to 1 


Rectangular Wing in Incompressible Flow 


The functions k-, ( s ) 


cir' 


k^Cs), and C ( k ) cir for rectangular wings 


of aspect ratio 4 and 6 have been calculated approximately and are tabulated 
in reference 7. The noncirculatory components of k-^s) are shown in refer- 


ence l8 to be 


0.860rt 

5T 


6 ( s ) and 


0 . 895 * 


&(s) for A = 4 and A = 6, 


respectively. In order to obtain the function cp(k), an exponential 
approximation was made for the function ^(s). These expressions are 
now presented, together with the expression for the vanishing-aspect-ratio 
case obtained from reference l6 and the corresponding oscillatory 
functions cp(k): 


M s )a_^ =1 (s) 


A=0 


(28) 


k 2 (s) A= 4 « 1.0 - 0.391e"°' 285s - 0.609e _:i - 658s 


(29) 


k 2 (s) A=6 « 1.0 - 0.535e'°* 2 " s - 0.465e- 2 - 00s 


(30) 
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9(k)^_o - 1.0 


(31) 


9(k) A= i, 


X _ 0-391 

,ik 0.285 + ik 


0.609 \ 

1.638 + ik/ 


(32) 




0-333 

0.299 + ik 


0.465 
2.00 + ik 


(33) 


Again, the asymptotic expressions could not be obtained except for the 
vanishing-aspect-ratio case: 


C(k) s 1 + 2ik 
|cp(k)| 2 = 1 


(k » 1.0) (34) 


Two-Dimensional Wing in Subsonic Compressible Flow 

For the two-dimensional wing in subsonic compressible flow, the 
following equations have been derived from references 2 and 3: 


k, (s) « 1.0 - 0.352e" 0,075 ^ s - 0.26le"°’ 5T2s + 0.669e" 1-89 ° s (35) 

-L M=0 . 5 

k-i (s ) , » 1.0 - 0.362e -0 - 06l4 ' 6s - 0.504e"°'^ 8ls + 0.7l4e'°- 958s (36) 

1 M=0 . 6 


k, (s) « 1.0 - 0.364e _0 -° 556s - 0.405e"°' 5 57s _ 0.4l9e“°- 9 ° 2s (37) 

1 M=0 . 7 

k 2 (s) M Q » 1.0 - 0.390e" 0 *° 7l6s - 0.407e -0 ’^ 7 ^ s - 0.203e -2 ' l69s (38) 

k 0 (s) . « 1.0 - 0.328e -O * O5l+5s - 0.430e _0 - 257s - 0.242e~ 1 - 1< - 6ls (39) 

2 V M=0.6 

k 2 ( s ) M= o * 1.0 - 0.402e _ °' 05 ^ 2s ~0.46le"°* 5125s - 0 . 137e _1 ^ 74s (40) 

cp(k) ~ ±t( 1 °- 39 ° - °- k01 ) 

; M=0.5 \ik 0.0716 + ik 0.37^ + ik 2.165 + ik/ 


( 41 ) 
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^^=0.6 ~ lk (ik “ 0. 


0.328 


0.430 


0545 + Ik 0.257 + ik 


“ lk (i - O.O^g 0 ! ^ 


0.461 


0.3125 + ik 


0.242 
1.461 + ik, 


_ 0 - 1-37 ) 

1.474 + ik/ 


(42) 


( 43 ) 


The function C(k) has been calculated from the coefficients com- 
piled in table I of reference 5* 

The asymptotic expressions can be shown to be 


C(k) ~ ^(l + i - - M 
v ' 2kM 

|cp(k)| 2 « -f 


rt^Mk 2 


(k » 1.0) 


( 44 ) 


provided it is assumed that the functions k^(s) and k 2 (s) and their 
first derivatives are continuous . ( See appendix B . ) 

As shown in appendix B, the determination of the asymptotic behavior 
is dependent not only on the continuity of the function and its deriva- 
tives, but also on the value of the function and its derivatives at 
s = 0. Therefore, although the k-^ and k 2 functions have been deter- 

2M 

mined numerically for s > - — — — , the known exact expression for 


s < 


2M 


1 + M 

as given in reference 19 may be utilized to obtain the 


1 + M 

asymptotic expressions for C(k) and cp(k), provided it is assumed that 
the functions k^s) and k 2 (s) and their first derivatives are 

continuous . 


Vanishing-Aspect-Ratio Delta Wing in Incompressible 
and Compressible Flow 

For a delta wing of vanishing aspect ratio in incompressible flow, 
the indicial lift functions have been obtained from reference l6 and are 


k L (s) =1 (s) +| 5 (s ) 

k 2 (s) 

k 2 (s) = 1.0 


(0 ^ s < 2) 
(s > 2) 


( 45 ) 

( 46 ) 
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The corresponding oscillatory functions are 


C(k) =1.0+| ik 

3 


( 47 ) 



2k 


(48) 


and 


C(k) » 1.0 + | ik 
3 


> 


(k » 1.0) (49) 



The k-^ function for the vanishing-aspect-ratio delta wing in com- 
pressible flow has been presented in reference 19* However, because of 
the difficulty of obtaining the oscillatory lift function C(k) from 
this function, the results derived in reference 10 for the oscillatory 
lift function for the delta wing are presented instead. The kg function 

has been derived in appendix C of this report and the corresponding cp(k) 
function was subsequently derived. The functions |cp(k)j^ and C(k) are 
presented in figures 9 and 10 for Mm = M ^ = 0.1 where m is the tangen - 

of the semiapex angle of the delta wing. The asymptotic behavior of 

|cp(k)|^ for Mm = 0.1 can be obtained if it is assumed that the kg(s) 

function and its first two derivatives are continuous. Based on the 
analysis obtained in appendix B and the kg function presented in appen- 
dix C, it can be shown that 


The asymptotic expression for the function C(k), also based on the 
assumptions stated previously, can be shown to be 



k 


(k » 1.0) 


(50) 


C(k) ~ 1 + 


jtmM mMk 


(k » 1.0) 


(51) 
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Two-Dimensional Wing in Sonic Flow 

The functions C I ^k ] _(s ) and C I ^k 2 (s) for a wing in two-dimensional 

sonic flow were obtained from the functions presented in reference 20 for 
a wing in two-dimensional supersonic flow by taking the limit as M — >1.0. 
The resulting expressions are 




\M s) = * 


2 \/s - 1 + cos ~ 1 ( p f ~ ) 


(0 ^ s ^ 1.0) 


(s > 1.0) 


(52) 


C I u k 2^ s ^ = 2s 


(0 < s ^ 1.0) 


'Ice 


kg(s) = | s cos-^S-g-S) + - (s > 1-0) 


(53) 


The function C L ^C(k) was obtained from reference 4 and Cj^q>(k) was 
derived by the use of equation (7). These expressions are 


Ct C(k) = e" lk + 4(1 + i)f(k) 

^ y2itk 


(54) 


C^k) = 2(1 - — - f(k) 
where the Fresnel integral f(k) is defined by 


f(k) = 


'k _ 


ix 


'0 )/2*x 


dx 


(55) 


(56) 


The unnormalized functions have been presented, since the theoretical 

value of Ct for two-dimensional sonic flow is infinite. 

% 


Two-Dimensional Wing in Supersonic Flow 

For a two-dimensional wing in supersonic flow, the following equa- 
tions have been derived from references 20 and 5 : 


NACA TN 3714-8 


15 


Ms) = I 


M 


0 < s < 


2M 


~ 1 + M 


k i (s> = £ 


1/rt , . _i 2M - Ms\ , 1 -1 2m + 

— - + sin - 1 - + - cos 

M\2 s / p 

R 


2 1 s - sM 2 


2M 


1 “ 1 / 


/-2M < s < -2M_) 

\M + 1 = M - 1/ 


k-, ( s ) = 1.0 


s > 


2M \ 

4 - 1 / 


k 2 (s) 

k 2 (s) 


2M 


M 


o < s < 

- “ M + 1 


( 57 ) 


-2- sfS + sin" 1 g£ - ^ + 1 cos" 1 2M 2 + s , sM 2 

2rtM \ 2 s / jt 2M 


2M < _2M_ 

M + 1 = M - 1 


( 58 ) 


k 2 (s) 


= 1.0 


s > 


2M 


M - 1 


C(k) = f 0 (M,o)) + 2ik[f 0 (M,a>) - f L (M,S)J 


where 


9(k) = f 0 (M,aj) 


2M k 

CD = — 


( 59 ) 

(60) 

(61) 


P 


and and are the Schwarz functions of order 0 and 1 , the Schwarz 

function of order n being defined as 


f n (M,S) = f 1 x"e-^J 0 (!f)ax 


(62) 


The functions f n (M,<jo) are tabulated for M > 1.0 in reference 21, 
for example. The asymptotic expressions are 


16 


NACA TN 3748 


CM - 1 


|cp(k)| 2 =* -K 


(k » 1.0) 


(63) 


4M 2 k 2 

Wide Delta Wing in Supersonic Flow 

For a wide delta wing in supersonic flow, the following equations 
have been derived from references 22 and 8: 


kl(s) ' “l 1 + 5?, 


= § 




w 


m(I ■ Wfe " ( k ' t) 


( 0< 

s < 

2M 


M + 1, 

cos -1 

2M 2 

- sM 2 


2M 

2M 

< s 

< 2M 


M 


k-^(s) = 1.0 


s > 


2M 


M - 1, 


(64) 


■ fc ° 2 




. / \ 1 ___-l 2M 2 - sM 2 + s P 2 ro „-l Ms - 2M _ 

= n cos 2M + 5*M s cos I 


§(§• “ !•) cosh -1 . ■ - s 

«\2 1 M|s - 2 


2M 


M + 


_ < s < 2M \ 

1 = M - 1/ 

M - V 


(65) 


k2(s) = 1.0 

C(k) = 2 | f 0 (M,oj) - f ^(M,aj) + ikjf 0 (M,aj) - 2f L (M,a3) + f £ (M,S)] j (66) 

(67) 


VM = i 


i ' 2ikf ° s#) ’ f ° (M,i5) 


3W 
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and 




!<P(k) 


1 2 _ p 2 (M 2 + l) 


l4M K 


(k » 1.0) 


(68) 


The function fJ — was calculated by numerical integration of equa- 
U \ M M / 

tion (62), inasmuch as no tabulations were available of the Schwarz 
function for M < 1. 


Wide Rectangular Wing in Supersonic Flow 

For a wide rectangular wing in supersonic flow, the following equa- 
tions have been derived from references 23 and 9: 


4 _ 2 ' 

v P 


M s) = M ‘ 


4 . 2 s - S 

M ^V S TT 


0 < s ^ 


2M 


M + 1 


(3 ' ^l) kl(s) 


4 

it 

1 

M 
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APPENDIX B 


ASYMPTOTIC BEHAVIOR OF OSCILLATORY LIFT COEFFICIENTS AS 
DETERMINED FROM INDICIA! LIFT FUNCTIONS 

In the following analysis it is shown that the asymptotic 
"behavior of the oscillating lift coefficients can he determined from the 
initial behavior and the discontinuities in the derivatives of the 
indicial lift functions . 

If \|r(k) represents either C(k) or cp(k) and if K(s) repre- 
sents the corresponding function k-j_(s) or k^Cs), then the reciprocal 

relation is 


(2) There is a sequence of points s.* (where j = 1, 2, . . .) 

at which one or more of the derivatives K' ' (s) (where 
n > N) has a finite discontinuity 

then by N + 1 successive integrations by parts, equation (7*0 can be 
expressed as 


For further integrations by parts, the discontinuities contribute terms 
of the form 



0 


( 7*0 


If K(s) has the following properties: 

(l) K(s) and all its derivatives up to and including K^(s) are 
continuous 




(76) 
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APPENDIX C 

INDICIAL LIFT FUNCTION FOR DELTA WING OF VANISHING 
ASPECT RATIO IN COMPRESSIBLE FLOW PENETRATING 
A SHARP -EDGED GUST 


The differential equation which governs the flow field for a very 
narrow delta wing has been shown in reference 19 to be 


0yy + ^zz 0tt 


(78) 


where the x, y, and z axes are fixed on the wing and the time variable 
is t 1 = t/a. The loading coefficient can be shown to be 


AP _ 4/1 d 

q " VIM ^t 


0 , 


(79) 


and the boundary conditions associated with this wing penetrating a 
sharp -edged gust are 


W u = -W 0 

w u = 0 


( x < Mt ) 

> 

(x > Mt) 


(80) 


where W u is the induced vertical velocity on the wing. As in refer- 
ence 19 the boundary conditions need only be satisfied over a span strip 
of the wing, since it was assumed that the velocity gradients in the 
y-, z-, and t-directions are independent of the gradient in the x-direction. 
See following sketch. 



Sketch (a) 
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Sketch (b) then represents the problem under consideration: 
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Sketch (b ) 


If the axes are transformed by the relation 


t = t 


x 

M 


(81) 
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then the problem can he represented by the following sketch: 



Sketch (c) 


In the lifting-surface analog, this corresponds to the problem of 
finding the velocity potential over a flat rectangular wing of low aspect 
ratio situated in a free stream at a Mach number equal to \[2. Solutions 
to this problem are given in reference 19 and are now presented; the sub- 
scripts represent the region under consideration. 

(0t)i = -W o ( 82) 
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emd, for large values of r (that is, for the higher number regions in 
sketch (c)), 


0T 


- -W Q f 



( 85 ) 


where an integral equation for 



this function is given in reference 19- 


together with a tabulation of 


Solving equations (82) to ( 85 ) for the velocity potential 0 and 
substituting into, equation ( 79 ) yields the following loading coeffi- 
cients for the various regions: 
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The corresponding lift coefficient per unit length is 
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Evaluation of equation (90) yields 



4Wcm T 
V b(x) 


(t < 2b(x) ) (91a) 





(t » 2b(x) ) (91b) 


where the subscripts indicate the regions in the xt-plane where Z(x) 

applies. See following sketch. 
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A plot 
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of the 


of equations (91a) 
seen that an error 

assumption made in 


and (91b) is presented in figure 21 where it 

of about 3 percent exists at — — — = 2 because 

b(x) 


equation (85). 


The lift coefficient for the wing is obtained from the equation 
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Substituting the appropriate expression from equations (91a) and (91h) 
into equation (92) gives 
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Evaluating equation (93) and making the substitution 
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Equation (94) has been evaluated numerically for Mm = ^ = 0>1 a™ 3 - the 

results are presented in figure 22, together with the results for M = 0 
as given in reference 17 - It might be of interest to note that for 
M = 0 equation (94) reduced to the equation given in reference 17 . 
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TABLE I— INDEX TO FIGURES, EQUATIONS, AND REFERENCES 


type of flow 

Wing plan form 

C(k) 

tp(k) 

C(k) = F(k) + IG(k) 
in figure - 

l |<p(k )| 2 

in figure - 

Incompressible 

Two dimensional 

Eq. (11) 

Eq. (12) 

1 

2 

Incompressible 

Elliptical 
A = 0 
A = 5 
A = 6 

A SS 00 

Eq. (21) 
Eq. (22) 

Eq. (23) 

Eq. (11) 

Eq. (24) 
Eq. (25) 
Eq. (26) 
Eq. (12) 

‘ ( a )3 

4 

Incompressible 

Rectangular 
A = 0 
A = 4 
A = 6 
A = oo 

> Ref. 16 
r Ref. 7 
Eq. (11) 

Eq. (31) 
Eq- (32) 
Eq. (33) 
Eq. (12) 

(a) 5 

6 

Subsonic compressible: 
M = 0 
M = 0.5 
M = 0.6 
M = 0.7 

Two dimensional 

Eq. (11) 
Uef. 3 

Eq. (12) 
Eq. (41) 
Eq. (42) 
Eq. (43) 

1 7 

- 8 

Compressible and 
incompressible : 
M = 0 

f = 0.1 

Vanishingly low 
aspect-ratio delta 

Eq. (4 7) 
Ref. 10 

Eq. (48) 
(b) 

) 9 

10 

Sonic 

Two dimensional 

Eq. (54) 

Eq. (55) 

11 

12 

Supersonic: 

M = f , f , 2, and f 

Two dimensional 

Eq. (59) 

Eq. (60) 

13 

14 

Supersonic : 

M - f , f , a, and f 

Wide delta 

Eq. (66) 

Eq. (67) 

15 

16 

Supersonic : 

M = i&, a, and ^ 

M = Jfi 

7 

M = 2 
3 

Wide rectangular 

JA-1 

< A = 2 
A = 4 

■ 

A = 1, 2, 4, and oo 

j 

Eq. (71) . 

Eq. (72) 

17(a) 

17(b) 

17(c) 

19(a) 

19(b) 

19(c) 

18(a) 

18(b) 

18(c) 

20(a) 

20(b) 

20(c) 


^Only circulatory component of lift plotted. 

Calculated numerically from k 2 (s) function given in appendix C. 



Figure 1.- The functions C(k) for a wing in two-dimensional incom- 
pressible flow. 
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Figure 2.- The function |cp(k)|^ for a wing in two-dimensional incom- 

pressible flow. 
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Figure 3 .- The functions C(k) cir for elliptical wings of various aspect 

ratios in incompressible flow. 
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Q 

Figure 4.- The functions |cp(k)| for elliptical wings of various aspect 

ratios in incompressible flow. 
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Figure 5.- The functions C(k) cir for rectangular wings of various aspect 

ratios in incompressible flow. 
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Figure 6.- The functions |ep(k)|^ for rectangular wings of various aspect 

ratios in incompressible flow. 
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Figure 7 .- The functions C(k) for a wing in two-dimensional subsonic 

flow. 
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Figure 8.- The functions |cp(k)|^ for a wing in two-dimensional subsonic 

compressible flow. 
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Figure 9 .- The functions C(k) for delta wings of vanishingly small aspect 
ratios in incompressible and compressible flow. 
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Figure 10.- The functions |cp(k)| 2 for a delta wing of vanishingly small 
aspect ratio in incompressible and compressible flow. 
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Figure 11.- The function CLC(k) for a wing in two-dimensional sonic flow. 
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Figure 12.- The function |c L cp(k)| 2 for a wing in two-dimensional sonic 

flow. 
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Figure 13 •- The functions C(k) 


for a wing in two-dimensional supersonic 
flow. 
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cp(k)| 2 for a wing in two-dimensional super- 
sonic flow. 


Figure 14.- The functions 
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Figure 15*- The functions C(k) 


for a wide delta wing in supersonic flow. 
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Figure l6.- The functions |cp(k)|^ for a wide delta wing in supersonic 

flow. 
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(a) A = 1. 

Figure 17*- The effect of Mach number on the functions C(k) for a wide 

rectangular wing in supersonic flow. 
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(b) A = 2. 


Figure 17 
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(a) A = 1. 


Figure 18.- 


The effect of Mach number on the functions |cp(k)| 
rectangular wing in supersonic flow. 
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Reduced frequency, k 

(b) A = 2. 

Figure 18.- Continued. 



NACA TN 37^8 


51 


|*00 



(c) A = 4. 

Figure 18.- Concluded. 
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(a) M = 10/7. 

Figure 19»- The effect of aspect ratio on the functions C(k) for a 
rectangular wing in supersonic flow. 


MCA TN 37^8 




Reduced frequency, k 


(b) M = 2. 

Figure 19 •- Continued. 
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(c) M = 10/3. 
Figure 19 .- Concluded. 
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Figure 20.- 


(a) M = 10/7. 

The effect of aspect ratio on the functions U(k 
a rectangular wing in supersonic flow. 
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(b) M = 2. 


Figure 20 


Continued 
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(c) M = 10 / 3 . 
Figure 20.- Concluded 
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Figure 21.- Growth of lift on a spanwise strip of a narrow delta wing 

penetrating a sharp-edged gust. 
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Figure 22.- The indicial lift function l^Cs) for a delta wing of 

vanishingly small aspect ratio in incompressible and compressible 
flow. 
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